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Abstract 

The hmitations of conventional imaging techniques have hitherto pre- 
cluded a thorough and formal investigation of the complex morphology of 
the left ventricular (LV) endocardial surface and its relation to the sever- 
ity of Coronary Artery Disease (CAD). However, recent developments in 
high-resolution Multirow-Detector Computed Tomography (MDCT) scanner 
technology have enabled the imaging of the complex LV endocardial surface 
morphology in a single heart beat. Analysis of high-resolution Computed To- 
mography (CT) images from a 320-MDCT scanner allows for the study of the 
relationship between percent Diameter Stenosis (DS) of the major coronary 
arteries and localization of the cardiac segments affected by coronary arterial 
stenosis. In this paper a novel approach for the analysis of the non-rigid LV 
endocardial surface from MDCT images, using a combination of rigid body 
transformation-invariant shape descriptors and a more generalized isometry- 
invariant Bag- of- Features (BoF) descriptor, is proposed and implemented. 
The proposed approach is shown to be successful in identifying, localizing 
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and quantifying the incidence and extent of CAD and thus, is seen to have a 
potentially significant clinical impact. Specifically, the association between 
the incidence and extent of CAD, determined via the percent DS measure- 
ments of the major coronary arteries, and the alterations in the endocardial 
surface morphology is formally quantified. The results of the proposed ap- 
proach on 16 normal data sets, and 16 abnormal data sets exhibiting CAD 
with varying levels of severity, are presented. A multivariate regression test 
is employed to test the effectiveness of the proposed morphological analysis 
approach. Experiments performed on a strict leave-one-out basis are shown 
to exhibit a distinct pattern in terms of the correlation coefficient within the 
cardiac segments where the incidence of coronary arterial stenosis is localized. 
Keywords: 

Ventricular endocardial surface, cardiovascular CT, non-rigid shape 
analysis, Shape-Index, Bag-of-Features. 

1. Introduction 

Since Coronary Artery Disease (CAD) is a leading cause of morbid- 
ity and mortality worldwide [7], techniques that increase the effectiveness 
and/or lower the costs of diagnostic or prognostic procedures associated with 
CAD are expected to have a significant clinical impact. CAD is caused by 
atherosclerosis or accumulation of lipoprotein plaque in the coronary arteries 
that supply blood to the myocardium or cardiac muscle tissue. Atheroscle- 
rotic plaques lead to the progressive narrowing or stenosis of the coronary 
arteries, resulting in reduced blood flow and consequently, reduced oxygen 
supply to the myocardium, a condition termed as myocardial ischemia. If 



2 



untreated, myocardial ischemia may lead to irreversible necrosis of the my- 
ocardium wherein the healthy myocardium is increasingly replaced by scar 
tissue, thus compromising the cardiac function and resulting in congestive 
cardiac failure. Some vulnerable plaques may suddenly rupture resulting in 
coronary artery occlusion and cardiac arrest, potentially leading to sudden 
death. 

X-Ray Coronary Angiography (XRA) is an invasive technique that is a 
clinically accepted standard for assessment of vascular morphology and the 
extent of vessel stenosis due to artherosclerotic plaque deposition. However, 
a comprehensive assessment of CAD necessitates a study of both, vascular 
morphology and cardiovascular function. Conventional cardiovascular func- 
tional assessment is performed via a stress-induced perfusion test that uses 
Magnetic Resonance (MR) or radionuclide Myocardial Perfusion Imaging 
(MPI). Since vascular morphology and cardiovascular function are imaged 
using separate modalities, the time and cost associated with a comprehen- 
sive assessment of CAD and the potential health risk to the patient associated 
with higher radiation doses are all significantly increased. 

Computed Tomography Coronary Angiography (CTCA) is a non-invasive 
imaging technique that allows for robust quantification of vascular morphol- 
ogy and has the potential for characterizing the atherosclerotic plaque com- 
position [21]. When performed using a 320-Multirow-Detector Computed 
Tomography (MDCT) scanner, a CTCA scan can be performed within a 
single heart beat, yielding images with an isotropic spatial resolution of 0.5 
mm in a volumetric fashion. The resulting images, in addition to providing 
vascular morphology information, are capable of providing significant details 
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Figure 1: Left ventricular trabeculae carneae (pointed by the red arrow) and papillary 
muscles (pointed by the yellow arrow) can be clearly seen in the CTCA images acquired 
using a 320-MDCT scanner. 

about the left ventricular (LV) endocardial surface structure. The CTCA 
images, backed by anatomical studies, have revealed that, instead of a sim- 
ple and smooth surface, the LV endocardial surface is composed of a complex 
structure of trabeculae carneae, which are small muscular columns that arise 
naturally from the inner surface of the ventricles (Figure [T]). Also revealed in 
the CTCA images are the papillary muscles which are attached to the cusps 
of the atrioventricular (i.e., the mitral and tricuspid) valves (Figure [T]). 

Alterations in the LV trabeculation structure have been clinically ob- 
served to closely associate with some cardiovascular diseases, such as my- 
ocardial noncompaction disease [1] and hypertrophy [7]. However, due to 
the limitations in the spatial resolution of conventional in vivo imaging tech- 
niques, very few research studies have been undertaken to formally investi- 



gate the LV trabeculation structure at a detailed level, and formally quantify 
the relationships between structural changes in LV trabeculation, changes in 
cardiovascular function and cardiovascular pathologies. If the analysis of 
the LV endocardial surface structure can be shown to provide significant in- 
sights into vascular morphology, cardiovascular function and the progression 
of CAD then it has the potential to provide significant additional or supple- 
mentary diagnostic value to the results of MPI, CTCA and XRA. In some 
cases, it could potentially reduce the need for invasive and/or stress-based 
testing procedures that are expensive, time consuming and pose a greater 
health risk to the patient, thus making CTCA imaging using MDCT scan- 
ners, a potential one-stop-shop technique for assessment of both, vascular 
morphology and cardiovascular function. 

The sub-endocardial layer of the myocardium is the first and the most 
susceptible cardiac region in the development of CAD on account of its in- 
herently higher oxygen consumption requirement and the restricted collat- 
eral blood flow resulting in reduced oxygen delivery during the progression 
of CAD. The presence of intermediate to severe CAD in patients is often 
clinically suggestive of varying degrees of myocardial ischemia in the cor- 
responding regions of the myocardium. Furthermore, the affected ischemic 
myocardium loses contractibility and has a tendency to get stiffer and be 
pushed outward by the high LV blood pressure. These myocardial changes 
are reflected in the LV endocardial surface morphology, in that the structure 
of the trabeculae and the papillary muscles on the LV endocardial surface are 
clinically observed to be substantially altered [20]. This clinically hypoth- 
esized association between the LV endocardial surface morphology and the 
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Figure 2: Illustration of the sequence of steps in the morphological analysis of the LV 
endocardial surface: (a) accurate mesh segmentation followed by (b) generation of a 17- 
segment LV surface model with demarcation of coronary arterial territories (red: LAD, 
green: LCX, blue: RCA), (c) feature vector generation and, (d, e) final generation of the 
BoF histogram via vector quantization (fc-means clustering). 
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incidence and severity of CAD constitutes the central thesis of this paper. 

The objective of this paper is to formally establish the clinically observed 
relationship between the incidence and severity of CAD, as quantified by the 
extent and location(s) of coronary arterial stenosis, and the structural alter- 
ations in the LV endocardial surface as observed in high-resolution MDCT 
images. The LV endocardial surface structure, specifically the structure of 
the LV endocardial trabeculae and papillary muscles, is formally character- 
ized using 3D morphological features obtained via 3D shape analysis algo- 
rithms. To the best of our knowledge, this research is amongst the first 
to formally characterize the LV endocardial surface structure and quantify 
the clinically observed association between the alterations in the LV endo- 
cardial surface structure and the extent and location(s) of coronary arterial 
stenosis. The research described in this paper could potentially add signif- 
icant value to the accuracy of diagnosis and effective management of CAD. 
It constitutes a first step towards formal investigation, quantification and 
verification of the clinically observed association between the LV endocardial 
surface morphology, cardiovascular function and vascular morphology which, 
in the long term, would aid in comprehensive assessment and understanding 
of the physiological mechanisms underlying the development and progression 
of CAD. 

2. Related Work 

Our preliminary work [20] demonstrated the potential diagnostic value 
of characterization of the LV endocardial surface structure in assessing the 
incidence of CAD. Although it produced some significant and encouraging 
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results, our preliminary work did exhibit shortcomings that stemmed pri- 
marily from the selection of shape descriptors used to characterize the LV 
endocardial surface. The two shape descriptors that were employed were 
based on the implicit assumption of rigidity of the LV endocardial surface 
as observed in the MDCT images. Consequently, the two shape descriptors 
were observed to have significant limitations in terms of their classification 
accuracy pO] . Although the MDCT image data were collected at a relatively 
steady phase of 75% in the R-R cardiac cycle, the continuous LV motion and 
the human error involved in the data collection, severely tested the limits of 
the rigidity assumption. It was clear that a more robust shape descriptor was 
needed to characterize the LV endocardial surface, i.e., one that is invariant 
to isometric global shape deformation [T5| 128]. 

Isometry-invariant shape descriptors have been the topic of recent re- 
search, especially in the context of content-based retrieval in shape databases |23] , 
shape symmetry detection [IB], dense surface correspondence determina- 
tion [28j and surface registration [24J. However, the problem in directly 
using existing shape descriptors developed for other applications, such as 
content-based image retrieval is that, these shape descriptors are designed 
primarily for differentiating between two distinct classes of objects (e.g., hu- 
mans in different poses versus dogs in different poses). In contrast, our goal 
is to differentiate objects within the same class, e.g., to classify whether a 
particular LV segment is normal or diseased, based on its surface geometry 
and surface morphology. 

Deriving a feature-based representation of the contents of an image gen- 
erally entails two stages of processing, feature detection and feature descrip- 
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tion [23]. The goal of feature detection is to localize relatively stable points 
or regions within an image that possess significant information and can be 
repeatedly and reliably detected in transformed versions of the image. Dif- 
ferent approaches are employed for feature detection based on the intrinsic 
nature of the features, i.e., point-based or region-based, and the desired scale 
of abstraction for the features. For example, in the case of the Scale Invari- 
ant Feature Transform (SIFT) feature detection technique, feature points of 
interest are localized by determining the local maxima of the discrete image 
Laplacian at varying scales of the Gaussian smoothing filter [10]. Point-based 
feature selection techniques could either select all the points in the image or 
perform a dense sub-sampling of the set of feature points [12] . 

Region-based feature detection techniques typically rely on segmentation 
of the cardiac images using shape priors [6l |25] . The high-dimensional shape 
prior is projected into a low-dimensional subspace using dimensionality re- 
duction techniques, e.g., variants of Principal Component Analysis (PCA), 
to constrain the shape variation. The detected features are tracked in a car- 
diac image sequence using a Kalman filter [6] or particle filter p5] whereas 
the shape prior is adapted using a learning algorithm such as AdaBoost [25] . 
Other region-based feature detection techniques for cardiac images include 
model-based deformable templates [5], Markov random fields [16], optical 
flow techniques [12] and combinations of the above [2]. 

The goal of the feature description stage is to arrive at a representation 
of the local image information in the neighborhood of the detected feature. 
For example, the SIFT technique represents each feature point by a 128- 
dimensional descriptor vector that is invariant to scaling. A similar approach 
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is taken in the Speeded- Up Robust Feature (SURF) transform approach [3], 
where a 64-dimensional descriptor for each feature point is computed effi- 
ciently using integral images. In order to achieve a parsimonious represen- 
tation, a feature vocabulary consisting of visual words is constructed by per- 
forming vector quantization in the feature descriptor space using a clustering 
technique such as the fc-means clustering algorithm. After having performed 
vector quantization, the individual feature descriptors are replaced by indices 
in the vocabulary of visual words. The aggregation of feature descriptors to 
describe the overall shape is achieved by generating a frequency histogram 
of the visual words in the vocabulary, termed as a Bag- of- Features (BoF). 

One of the prominent implementations of the BoF concept is Video Google, 
a web-based application designed by Sivic and Zisserman [26], ET] for object- 
based search in large image and video collections. In Video Google an image 
is represented by a collection of visual words. The visual words are generated 
using feature detectors and feature descriptors and indexed in a visual vo- 
cabulary. The BoF representation is constructed from the frequencies of the 
visual word occurrences in the image. Images with similar visual informa- 
tion are observed to have similar BoF representations. Thus, the comparison 
of BoF representations allows one to retrieve similar images. Likewise, two 
shapes are compared by comparing their respective visual word frequency his- 
tograms or BoF representations. This reduces the shape similarity problem to 
the problem of frequency histogram comparison. Shape Google, a technique 
for shape-based search in large image collections developed by Ovsjanikov et. 
al. [23] and, the shape comparison approach proposed by Toldo et. al. [29] 
are amongst the prominent works on BoF-based shape description and shape 
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comparison, and the work described in this paper is influenced by both of 
these papers [23], [29] . 

As an extension to the general BoF approach, we propose a novel BoF- 
based shape analysis approach designed specifically for cardiac image anal- 
ysis, and for medical image analysis in general. The proposed BoF-based 
shape analysis approach is designed to quantify the relationship between the 
incidence, severity and localization of CAD and the structural alterations in 
the LV endocardial surface. While several works have proposed feature-based 
approaches for characterization of rigid 3D shapes, very few are capable of 
dealing with non-rigid 3D shape deformations [T^ . To the best of our knowl- 
edge, this is one of the first works that uses a BoF-based shape analysis 
approach for comparison of non-rigid deformable 3D shapes in the context 
of cardiovascular imaging in particular, and medical imaging in general. 

In this paper, a BoF-based approach is proposed to encapsulate the local 
and global geometry, local surface orientation and global contextual infor- 
mation for the LV endocardial surface. The proposed approach is shown 
to result in a robust feature vector for the purpose of morphological anal- 
ysis of the LV endocardial surface. As shown in the experimental results, 
the proposed approach results in successful localization of coronary arterial 
stenosis which serves to strengthen the clinically observed relationship be- 
tween the incidence and severity of CAD and morphological alterations in the 
LV endocardial surface. The sequence of steps in the proposed approach for 
morphological analysis of the LV endocardial surface is depicted in Figure |2] 

The remainder of the paper is organized as follows. In Section [3} the 
main contributions of the paper are outlined; in Section |4| the proposed LV 
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surface segmentation and LV shape analysis procedures are detailed; and in 
Section [5| experimental results are presented. Finally, in Section [6| the paper 
is concluded with a brief discussion of the proposed approach and an outline 
of directions for future work. 

3. Contributions of the Paper 

The paper makes two primary contributions as described below: 

1. The Bag- of- Features (BoF) framework for non-rigid shape analysis is 
adapted for the purpose of cardiac shape analysis which is an impor- 
tant problem in cardiac imaging in particular and medical imaging in 
general. It is also important to note that the proposed approach is 
sufficiently general to be applicable to other problems in medical imag- 
ing such as 3D shape analysis of the human brain and 3D shape-based 
search and retrieval in large medical databases. 

2. The paper proposes a geometric and machine learning-based model of 
the relationship between localized changes in the LV endocardial sur- 
face morphology and the incidence and extent of stenosis in specific 
coronary arteries. This constitutes an important initial step towards 
clinical understanding of the complex relationship between coronary 
arterial stenosis and its effect on the morphology of the LV endocar- 
dial surface. To the best of our knowledge, this paper represents one 
of the first attempts to model this complex clinical relationship in a 
mathematically structured manner. 



12 



4. MDCT Image Segmentation and LV Shape Analysis 

4.I. LV endocardial surface segmentation and meshing 

The trabeculae carneae along the LV endocardial surface can be broadly 
classified into three different morphological types: (a) those that lie along 
the entire length of the LV wall forming prominent ridges; (b) those that are 
fixed at their extremities but free in the middle; and (c) those that connect 
the roots of the papillary muscles to the ventricular wall. These different 
trabeculae morphologies result in a complex LV endocardial surface topology. 

A 3D level set approach is employed to segment the LV endocardial sur- 
face while adapting to the topological changes caused by the complex trabec- 
ulation structure. A median filter-based denoising procedure is employed on 
the 3D MDCT data prior to segmentation in order to suppress noise while re- 
taining the edges in the MDCT images. Unless mentioned otherwise, the size 
of the median filter is empirically set to 7 x 7 based on the MDCT data set. A 
level set-based segmentation procedure without reinitialization, as proposed 
by Li et al. [9], is applied to the median-filtered 3D image data set followed 
by the Marching Cubes procedure [11] to generate the surface meshes. The 
surface meshes are subsequently denoised via a mean face normal filtering 
procedure proposed by Zhang and Hamza to obtain the smooth shape 
of the LV myocardial surface. 

4-2. Data preparation 

The standardized myocardial segmentation model proposed by the Amer- 
ican Heart Association (AHA) is adopted to facilitate accurate segmentation, 
understanding and localization of cardiac anatomy and pathology [H]. The 
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Figure 3: Bull's eye view of the 17-segment AHA model. The numerical label for each LV 
segment is color coded to denote the territory of the major coronary artery it falls under, 
i.e., Left Anterior Descending (LAD): yellow; Right Coronary Artery (RCA): cyan; and 
Left Circumflex Artery (LCX): black. 

AHA-approved 17-segment cardiac model P] is adapted to divide the left 
ventricle into 17 segments for more accurately localized shape analysis (Fig- 
ure [3]). The long axis of the left ventricle is first computed to divide the 
LV endocardial surface into 4 main parts along the longitudinal orientation 
termed as the apex, apical, mid-cavity and basal. Further division of the LV 
endocardial surface along the short axis view is tackled by exploiting knowl- 
edge of cardiac anatomy. Three landmark points are considered across the 
septum to carry out the division process. This results in the division of apical 
segment into 4 parts and the mid- cavity and basal segments each into 6 parts 
whereas the apex remains undivided. At the end of the segmentation process, 
the LV endocardial surface is divided into 17 segments in conformity to the 
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AHA model as shown in Figure [3] 

4-3. LV endocardial surface descriptors based on the rigidity assumption 

In our preliminary work [20], we considered two primary shape descrip- 
tors, to characterize the LV endocardial surface morphology, i.e., the D2 
shape descriptor, a global shape descriptor proposed by Osada et al. |22] and 
the shape index, a local shape descriptor first introduced by Koenderink |8] 
and subsequently modified by Zaharia and Preteux [30j. Both descriptors 
make an implicit assumption of rigidity of the endocardial surface and are 
invariant to scale and 3D rigid body transformation (i.e., rotation and trans- 
lation) in Euclidean space. 

The D2 shape descriptor is a global shape signature of a 3D object and is 
represented by a probability distribution function. The probability distribu- 
tion function is obtained via sampling of a pre-specified local shape function 
which measures a basic geometric property of the underlying 3D shape. In 
the case of the D2 shape descriptor, the shape function is the distance be- 
tween two randomly sampled vertices on the mesh-based representation of 
the 3D LV surface. Osada et al. [22] have shown that the simple D2 shape 
descriptor, when modeled as a probability distribution function over the en- 
tire 3D object surface, can serve as an effective object signature for a variety 
of applications such as 3D shape-based retrieval from a 3D object database. 

The shape index, in contrast, describes the local shape of the 3D surface 
mesh based on the values of the principal surface curvatures computed within 
a local neighborhood of a given 3D surface point. The shape index is partic- 
ularly effective in describing the local extrinsic geometry of the surface. The 
shape index Ip of a surface point p is a function of the two principal surface 
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curvatures Ki{p) and k,2{p) associated with point p and is defined as follows: 



where, ni{p) > i^2{p)- 

The value of the shape index lies within the interval [0, 1] and enables 
surface-based representation of basic elementary shapes such as convex, con- 
cave, rut, ridge and saddle [30j . Since the LV endocardial surface is observed 
to exhibit both, global structure as well as detailed local structure and since 
there have been no previous studies detailing which of the shape properties 
(local or global) convey more valuable clinical information, we considered 
both local and global shape descriptors in our preliminary work [20j. How- 
ever, the results of our preliminary work showed the shape index to be more 
effective than the D2 shape descriptor for classification of the LV endocardial 
surface samples into diseased and non-diseased classes. 

4-4- LV endocardial surface descriptors based on non-rigid deformation 

BoF feature-based approaches for shape representation typically consist 
of two steps, i.e., feature detection and feature description. For representa- 
tion of object shapes encountered in common imagery, the feature detectors 
and feature descriptors are designed to be invariant under affine transfor- 
mation to take into account the different possible viewpoints under which 
a 3D object may appear in an image. For non-rigid characterization of the 
LV endocardial surface, the type of invariance required is different from the 
one required when dealing with common images since the number of surface 
transformations in the case of the former is much larger, including bending, 




(1) 
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changes in pose and changes in connectivity. In our case, the LV endocardial 
surface deformation can be constrained to he within a pre-specified range 
since all the CTCA images are acquired at a relatively steady phase (75%) 
within the R-R cardiac cycle. The aforementioned steady phase assump- 
tion allows the LV surface deformation to be approximated by an isometric 
deformation. Consequently, the LV endocardial surface descriptors need to 
be robust enough to capture the important characteristics of the underlying 
shape while exhibiting invariance to isometric deformation. Since the LV 
endocardial surface is typically less rich in terms of features than common 
images, the task of detecting a large number of stable, repeatable and iso- 
metric deformation-invariant feature points is significantly harder. In this 
work, the features selected for description of the LV endocardial surface are 
inspired by the work of Toldo et. al. [29] . 

4.4- i- Feature detection 

Although feature detection is not mandatory, it can improve computa- 
tional efficiency in practice. In our case, we have chosen dense random sam- 
pling as the feature detection strategy. We considered as feature points, 500 
randomly sampled surface points from each of the 17 LV endocardial surface 
segments. The underlying motivation was to ensure that the sampling is as 
dense as possible so that the entire spectrum of interesting points on the LV 
endocardial surface is considered for further analysis. In situations where 
availability of computational power is not a constraint, the feature selection 
step can be ignored and all the points on the surface can be considered for 
the ensuing feature description step. 
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Figure 4: Illustration of the relationship between the (a) shape index Ip and (b) curvedness 
Cp surface descriptors for a sample LV endocardial surface. 

4 ■4- 2- Feature description 

The feature descriptors are designed to ensure that both the local and 
global geometric properties of the surface are represented in sufficient detail. 
Inspired by the work of Toldo et. al. [29j, four surface descriptors are con- 
sidered to represent each sampled surface point. The first three descriptors 
represent the local geometry of the surface around a specific point whereas 
the fourth is a contextual descriptor which encapsulates the global properties 
of a specific point with respect to other points on the surface. The descriptors 
are detailed below: 

4-4-S-i- Shape Index {Ip). The shape index Ip of a surface point p, is com- 
puted using (equation Q) as previously described. The shape index is partic- 
ularly important for its ability to quantitatively characterize the local shape 
of a surface. 
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4-4-3-2. Curvedness iCp). The curvedness descriptor represents the degree 
of local curvature of a surface. Normalized curvedness also encapsulates 
the rate of change of local surface curvature: a value close to 1 implies a 
very gentle change, whereas a value close to implies a very sharp change. 
Generally, points on a smooth spherical object have small curvedness values. 
The curvedness Cp at a surface point p, as proposed by Koenderink [8], is 
defined as a function of the two principal surface curvatures ni{p) and K2{p) 
associated with surface point p as follows: 



where, Ki{p) > K2{p)- The relationship between the two surface descriptors, 
i.e., the shape index Ip and the curvedness Cp for a typical LV endocar- 



curvedness Cp is also invariant to scale and 3D rigid-body transformation 
(i.e., translation and rotation) in Euclidean space. 

4 -4 ■2-3. Normal Orientation [Op). It has been clinically observed that the in- 
cidence of myocardial infarction and myocardial ischemia, can cause changes 
in the 3D orientation of the LV trabeculation structure due to the build up 
of scar tissue. Consequently, the local 3D orientation of the LV endocardial 
surface is incorporated within the feature vector. The normal orientation Op 
at surface point p is defined as the angle between the unit normal vector at 
surface point p and the XZ-plane. 

4-4-2-4- Geodesic Contextual Descriptor [GCD'^). The geodesic contextual 
descriptor is intended to provide a means for describing the overall shape 





dial surface is depicted in Figure [4j Like the shape index Jp, the value of 
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and measuring shape similarity based on point correspondences. Contextual 
information is used to describe a specific surface point p with respect to the 
3D shape as a whole. The geodesic contextual descriptor for a surface point p 
encapsulates the distribution of relative positions of other points on the same 
surface with respect to p, thus summarizing global shape in an informative 
and, most importantly, in an isometric deformation-invariant manner. 

Moreover, since BoF-based methods typically employ a frequency dis- 
tribution histogram, there is inherent loss of contextual information. An 
example of this is a Google search based on the keyword Matrix that yields 
mutually distinct results; results pertaining to the movie titled Matrix and 
results pertaining to the concept Matrix from linear algebra. These semanti- 
cally distinct results can be disambiguated only via incorporation of contex- 
tual information. Since it is difficult to incorporate contextual information 
within a frequency distribution histogram, it needs to be incorporated within 
the feature vector. The geodesic contextual descriptor, denoted by GCD^, 
is an isometric deformation-invariant contextual descriptor that encapsulates 
the relative position of a surface point p in relation to the other points on the 
surface segment. It is characterized by a 20-bin histogram which is generated 
by computing the normalized geodesic distance between the surface point p 
and other sampled points on the surface segment. The geodesic contextual 
descriptor GCD"^ is invariant to scale, 3D rigid body transformation and 
isometric deformation. 

The feature detection and feature description procedures described above 
result in a 23-tuple feature vector for each surface point p denoted by Fp = 
(/p, Cp, 9p, GCDf) as depicted in Figure |2| 
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Figure 5: Concept visualization for the BoF histogram generation procedure. The fre- 
quency distribution histogram of the sampled visual words is illustrated for a sample LV 
endocardial surface. 

4-4-3- Construction of visual vocabularies 

The feature vectors Fp at each surface point are clustered in order to 
obtain the visual words. Assuming that the local descriptors are computed 
for a set of stable surface points, the feature vector space is quantized to 
obtain a compact representation for the vocabulary of visual words, in a 
manner similar to the Shape Google approach [23]. A vocabulary of visual 
words is defined as a set of representative vectors in the descriptor space 
or feature space, obtained by means of unsupervised learning, i.e., vector 
quantization via fc-means clustering in our case. More formally, a vocabulary 
of visual words is defined as a collection V = {vi, . . . ,Vk} where Vi is the 
centroid of the ith cluster and the clusters represent the visual words. We 
have experimentally chosen /c = 20 in the fc-means clustering algorithm for 
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generating the final histogram and all the subsequent experiments are based 
on this particular value of k. Figure |5] depicts the frequency distribution 
histogram of the sampled visual words for a sample LV endocardial surface. 

5. Experimental Results 

5.1. Experimental setup 

The proposed methods for segmentation, meshing and endocardial surface 
shape description were tested on 32 MDCT data sets consisting of 16 data 
sets from cardiac patients and 16 data sets from normal subjects. The inci- 
dence of single- or multi- vessel obstructive CAD was found in the three major 
coronary arteries, i.e., the Left Anterior Descending Artery (LAD), Left Cir- 
cumflex Artery (LCX) and Right Coronary Artery (RCA) using XRA, and 
further confirmed by Myocardial Perfusion Imaging (MPl) and Fractional 
Flow Reserve (FFR) tests performed on the patients. The individual terri- 
tories for each of the three major coronary arteries are depicted in Figure [3] 
in conformity to the American Heart Association (AHA) convention. 

The cardiac patients and normal subjects were subject to a contrast- 
enhanced CTCA scan on a 320-MDCT scanner using a standard protocol 
with electrocardiogram (ECG) gating. The resulting images were recon- 
structed at a relatively steady state of 75% in the R-R cardiac cycle to ensure 
minimal ventricular motion. This ensures that the endocardial surfaces are 
reconstructed at a relatively fixed cardiac dilation stage, and the subsequent 
shape analysis is minimally affected by cardiac motion. This method still 
results in the presence of motion artifacts in the acquired images due to var- 
ious factors such as the high velocity of the ventricular motion and irregular 
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ECG. However, the assumption of isometric deformation enables us to bet- 
ter account for these motion artifacts. The segmentation method described 



in Section |4.1| was used to generate topologically correct and geometrically 
accurate data for subsequent analysis. 

The experimental results are detailed for three different scenarios. In 
the first scenario, the problem is modeled as one of global classification 
wherein the shape index (Jp) histogram, the rigid-body shape descriptor 
(Jp, Cp, 9p) histogram, the combined rigid-body and isometric deformation- 
invariant shape descriptor (Jp, Cp, dp, GCDp^) histogram and the BoF his- 
togram are each used to classify the entire LV endocardial surface into one of 
two classes, i.e., normal and diseased. In the second scenario, the BoF-based 
approach is used to model the problem as one of localized classification where 
each LV segment is labeled as normal or diseased based on the extent of di- 
ameter stenosis (DS) in the major coronary arteries supplying blood to the 
LV segment under consideration. The LV segment is considered normal if the 
percent DS value of the corresponding major coronary artery is less than 70% 
and diseased if it is > 70%. In the third scenario, the BoF-based approach is 
used to model the problem as one of multivariate regression where the exact 
percent DS values of each of the major coronary arteries are treated as the 
labels based on which the correlation coefficients for each of the LV segments 
are computed. 

5.2. Segmentation results 

The results of the LV endocardial surface segmentation have already been 
proven to be reasonably accurate in our preliminary study [20]. The spatial 
distribution of the trabeculation was observed by clinicians to vary with the 
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Figure 6: Illustration to compare the quality of the mesh generated by our segmentation 
algorithm. In particular, the three colors of the LV endocardial surface show the smooth- 
ness of the surface for different sizes of median filter kernel, i.e., 5x5 (purple), 7x7 (sky) 
and 9 X 9 (blue). 

location within the left ventricle; thus providing the rationale for using the 
standard 17-segment AHA model to perform localized shape analysis. Fur- 
thermore, the proposed segmentation approach has already shown a visually 
observable distinction in trabeculation between normal and diseased hearts, 
yielding a classification accuracy greater than 80% with simple rigid-body 
surface descriptors and a nearest-neighbor classifier [20]. Our previous work 
has demonstrated the accuracy of the proposed segmentation method as well 
as its applicability for subsequent quantitative shape analysis [20]. In this 
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particular experiment, we varied the size of median filter kernel to determine 
its effect on the final surface mesh. The variation of the smoothness of the 
LV endocardial surface with respect to three different median filter kernel 
sizes, i.e., 5 x 5, 7 x 7 and 9x9 are depicted in Figure |6j In particular, for all 
the later experiments, we fixed our median filter kernel size at 7 x 7 as it was 
observed to provide a good balance between smoothness and preservation of 
geometric details of the LV endocardial surface. 

5.3. Rigid shape description results 

The experimental results from our preliminary study [20j demonstrated 
the effectiveness and limitations of both, the D2 shape descriptor and shape 
index, for describing the LV endocardial surface. The results suggested that 
the information derived from the D2 shape descriptor is inadequate for distin- 
guishing between normal and diseased hearts. Since the D2 shape descriptor 
is a global shape descriptor designed specifically to distinguish between ob- 
jects of different classes, e.g., cup versus car, it does not capture adequately 
the local spatial details of the underlying shape. On the other hand, since 
the shape index Ip is a local shape descriptor, it successfully captures the 
local shape details. A 20-bin shape index histogram was generated for each 
LV segment. In particular, the difference in shape between the normal and 
diseased hearts was represented and visualized in 17 x 20 dimensions (17 
segments per left ventricle x 1 histogram per LV segment x 20 bins per 
histogram) . 

A typical binary classifier can be trained to separate the normal training 
samples from the diseased samples in high-dimensional feature space by de- 
termining a projection function that maximizes classification accuracy. Lin- 
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(a) (b) (c) 

Figure 7: Illustration to compare the results generated by the D2 shape descriptor (b) 
and the shape index (c) for the mid-cavity segment 11 (blue) of the left ventricle (a). 12 
sample datasets are considered here where the 6 normal left ventricles are represented in 
red and the 6 diseased ones are represented in blue in the histograms in (b) and (c). 

ear Discriminant Analysis (LDA) is a projection and classification method 
that maximizes the between-class scatter defined by X]j=i(/^i ~ ~ lA^ 
and minimizes the within-class scatter defined by X]j=i " ~ 

/ij)^ where /i is the mean of all classes, is the zth sample of the jth class, 
yUj is the mean of jth class, Nj is the number of samples in the jth class and 
c is the number of classes. The 17 x 20 = 340-dimensional feature vector 
generated from the Shape Index Jp-based histogram is projected to a sub- 
space of c — 1 dimensions using LDA. In our case c = 2 since the problem 
is one of classification of the samples into one of two classes, i.e., normal or 
diseased. Classification is done using a fc-nearest neighbor (A;-NN) scheme in 
the resulting ID subspace. 

We tested the LDA-based fc-NN classifier for k = 1 and 3, and obtained 
the same result: 26 out of 32 samples were classified correctly resulting in 
an overall accuracy of 81.25%. Table [T] shows the confusion matrix for the 
diagnosis accuracy. The false alarm rate (normal sample classified as dis- 
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eased) was observed to be 18.75% whereas the miss rate (diseased sample 
classified as normal) was observed to be 18.75%. The experimental results 
confirmed that for identification and localization of cardiovascular patholo- 
gies, local shape descriptors convey much more information than global ones. 
The binary classification experiment was repeated for the histograms gener- 
ated by the feature vectors (Jp, Cp, Op) and (Jp, Cp, Op, GCD"^) and the results 
tabulated in Tables |2] and [3} respectively. 

Table 1: Confusion matrix to illustrate the classification accuracy of the Shape Index 
/p-based histogram. 





Classified Diseased 


Classified Normal 


Actual Diseased 


13 


3 


Actual Normal 


3 


13 



Table 2: Confusion matrix to illustrate the classification accuracy of the Shape Index, 
Curvedness, Normal Orientation (/p, Cp, 0p)-based histogram. 





Classified Diseased 


Classified Normal 


Actual Diseased 


13 


3 


Actual Normal 


3 


13 



The Shape Index, Curvedness, Normal Orientation (/p, Cp, ^p)-based his- 
togram was observed not to result in any improvement in classification over 
the Shape Index Jp-based histogram (Tables [l] and In the case of the 
(Jp, Cp, Op, CCDp°)-based histogram, 27 out of 32 samples were classified cor- 
rectly resulting in an overall accuracy of 84.37%. Table[T]shows the confusion 
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matrix for the diagnosis accuracy. The false alarm rate (normal sample classi- 
fied as diseased) was observed to be 18.75% whereas the miss rate (diseased 
sample classified as normal) was observed to be 12.5%. Thus, the incor- 
poration of the 20-bin Geodesic Context Descriptor GCD^^) in the feature 
vector was observed to improve the classification accuracy on account of its 
invariance to isometric deformation. 

Table 3: Confusion matrix to illustrate the classification accuracy of the 23-tuple feature 
vector {Ip, Cp,9p, GCDp°)-hased histogram. 





Classified Diseased 


Classified Normal 


Actual Diseased 


14 


2 


Actual Normal 


3 


13 



To demonstrate the superiority of the BoF-based description, the bi- 
nary classification experiment described above was repeated for the BoF 
histogram. In this experiment, 29 out of 32 samples were classified correctly 
resulting in an overall accuracy of 90.62%. Table |4] shows the confusion ma- 
trix for the diagnosis accuracy. The false alarm rate (normal sample classified 
as diseased) was observed to be 6.25% whereas the miss rate (diseased sample 
classified as normal) was observed to be 12.5%. 

Table 4: Confusion matrix to illustrate the classification accuracy of the BoF histogram. 





Classified Diseased 


Classified Normal 


Actual Diseased 


14 


2 


Actual Normal 


1 


15 
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Bag-of-Words Histogram 




Visual Words 



Figure 8: Illustration of a sample BoF frequency histogram which demonstrates its dis- 
criminative power. The normal (blue) subjects show a very distinct frequency histogram 
pattern when compared to the diseased subjects (green). A particular area of the BoF 
frequency histogram is shown here to clarify the dominant presence of diseased subjects 
over normal subjects in this particular region which, in turn, is reflected in the higher 
classification accuracy for the BoF histogram. 

5.4- Non-rigid shape descriptor-based localized classification results 

A series of experiments were performed to demonstrate the effectiveness 
and limitations of the non-rigid shape descriptor Fp = {Ip,Cp,9p,GCDp^). 
In this series of experiments, a coronary artery was considered as diseased 
or stenotic if the extent of stenosis is 70% or greater. The LV myocardial 
segments were labeled as diseased by a cardiologist if they were supplied 
blood by stenotic arteries. The already available percent Diameter Stenosis 
(DS) data were used for determining whether a coronary artery is normal or 
diseased. An artificial neural network (ANN), employing a multilayer per- 
ceptron (MLP) architecture with a single hidden layer and a learning rate of 
0.3, was used for the purpose of classification of the LV segments in a manner 
similar to that described in [2T] . The 20-bin BoF frequency histograms, gen- 
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Figure 9: Illustration of the classification accuracy for detection of coronary arterial steno- 
sis based on the change of surface morphology in the 17 LV segments in the AHA model. 
Higher gray values denotes higher classification accuracy and vice versa. The numerical 
label for each LV segment is color coded to denote the territory of the major coronary 
artery it falls under: LAD: red; RCA: blue; and LCX: green. 



erated via the vector quantization procedure (Section 4.4), for a particular 
LV segment from all the LV datasets were used as the inputs to the MLP 
ANN. Figure |8] illustrates the discriminative power of the BoF frequency his- 
tograms in being able to distinguish between normal and diseased subjects. 
The ANN-based classification procedure was carried out within a strictly 
leave-one-out setting. The output of the MLP ANN is whether a particular 
LV segment can be classified as normal or diseased. 

The success rate for detection of stenosis for a specific coronary artery 
is shown in Figure [9j The classification results depict a clinically observed 
relationship between the coronary arterial stenosis and the affected segment 
in the 17-segment AHA model. The lower classification accuracy in the 
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basal area (segments 1-6) can be explained by the clinical observation that 
several instances of coronary arterial stenosis encountered in this study are 
located in the mid to distal portion of the coronary arteries that only effect 
the mid-cavity (segments 7-12) and apical (segments 13-16) portions of the 
LV endocardial surface. Furthermore, another probable reason for the lower 
classification rate in the basal area is that the apical and mid-cavity segments 
exhibit greater endocardial trabeculation structure than the basal segments, 
which translates to more reliable endocardial surface morphology information 
that can be used for the purpose of classification in the case of the apical and 
mid-cavity segments. 

5.5. Non-rigid shape descriptor-based localized multivariate regression results 



In this experiment, the classification problem of Section 5.4 was recast as 
a multivariate linear regression problem wherein the exact percent DS data 
was considered as the ground truth. The regression procedure was carried 
out within a strictly leave-one-out setting. Multivariate linear regression 
attempts to fit a linear model to multiple (i.e., more than one) independent 
(or explanatory) variables xi,X2, ■ ■ ■ Xn, to obtain an estimate for a single 
dependent (or response) variable y as shown below: 

y = l3o + (3iXi + P2X2 ■■■ + PnXn + e (3) 

where the /3j's are the model coefficients and e is the fitting error. In our 
case, the dependent variable y denotes the exact percent DS value and the 
independent variables Xj, i = 1, . . . ,n = 20 denote the components of the 
20-tuple BoF histogram feature vector generated via the vector quantization 



procedure described in Section AA During the training process, the exact 
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percent DS data for a major coronary artery (LAD, LCX and RCA) and the 
20-tuple BoF histogram feature vectors corresponding to the LV segments 
that comprise the territory of that coronary artery are fed as inputs to the 
multivariate hnear regression procedure. The multivariate linear regression 
procedure determines the model coefficients, i.e., /3j's, that minimize the 
mean squared error given by ^S^^^e^ where m is the number of exact percent 
DS value points associated with the major coronary artery (equation (|3])). 

During the testing procedure, the estimated DS percent value, i/est, for 
the major coronary is computed using the 20-tuple BoF histogram feature 
vector (a;i, X2, . . . , X20) and the coefficients (/3o, /^i, . . . , f32o) as follows: 

Vest = /3o + Pi^i + 132X2 ■ ■ ■ + l3 nXn (4) 

The correlation coefficient is computed between the estimated outcome Uest 
of equation ^ and the exact percent DS data y (i.e., the ground truth). 

The correlation coefficient for detection of morphological changes (as 
quantified using the 20-tuple BoF histogram feature vectors) with respect 



to the percent DS values in a localized LV segment is shown in Figure 10 
The regression results were seen to depict an intriguing and clinically ob- 
served relationship between the coronary arterial stenosis and the affected 
LV segment in the 17-segment AHA model. The lower correlation coefficient 
values in the basal area (segments 1-6) compared to those in the mid-cavity 
(segments 7-12) and apical portions (segments 13-16) of the LV endocardial 



surface can be explained in a manner similar to that described in Section 5.4 
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Figure 10: Illustration of the correlation coefficient for detection of coronary artery stenosis 
based on the change in surface morphology in the 17 LV segments in the AHA model. A 
bluish color denotes a higher correlation coefficient value whereas a reddish color denotes 
a lower correlation coefficient value. The numerical label for each LV segment is color 
coded to denote the territory of the major coronary artery it falls under: LAD: yellow; 
RCA: green; and LCX: black. 

6. Discussion and Conclusion 

To the best of our knowledge, this paper is one of the earliest works that 
studies the relationship between coronary artery stenosis and the morpholog- 
ical alterations in the LV endocardial surface using high-resolution MDCT 
data, and demonstrates its potential predictive value for the diagnosis of the 
incidence and severity of CAD. One limitation of our work is the use of the 
LV territories defined on the basis of the current American Heart Associa- 
tion (AHA) convention. However, the AHA convention may not accurately 
reflect the true anatomy in individual patient cases which in turn can be an 
explanation for the low classification accuracy in the localized classification 
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results (Section 5.4). 



Our investigation also sheds new light on the localization of LV regions 
that are the most affected by coronary artery stenosis, a phenomenon which is 
yet to be fully explained. The association between the morphological features 
of the LV endocardial surface and cardiovascular function will be explored 
in our future work. In particular, we aim to investigate the correlation be- 
tween the LV endocardial surface morphology and the results of Myocardial 
Perfusion Imaging (MPI) and Fractional Flow Reserve (FFR) tests in addi- 
tion to the coronary arterial stenosis results obtained via X-Ray Coronary 
Angiography (XRA). 

References 

[1] Y. Agmon, H.M. Connoll, L.J. Olson, B.K. Khandheria, and J.B. Se- 
ward, Noncompaction of the ventricular myocardium. Jour. American 
Soc. Echocardiography, Vol. 12(10), Oct. 1999, pp. 859-863. 

[2] K. Althoff, G. Hamarneh, and T. Gustavsson, Tracking contrast in 
echocardiography by a combined snake and optical flow technique, IEEE 
Computers in Cardiology, Vol. 3653, Cambridge, MA, Sept. 2000, pp. 
29-32. 

[3] H. Bay, A. Ess, T. Tuytelaars, L. Van Cool SURF: speeded up ro- 
bust features. Computer Vision and Image Understanding (CVIU), Vol. 
110(3), 2008, pp. 346-359. 

[4] M. D. Cerqueira, N. J. Weissman, V. Dilsizian, A. K. Jacobs, S. Kaul, 
W. K. Laskey, et al. Standardized myocardial segmentation and nomen- 

34 



clature for tomographic imaging of the heart, Circulation, Vol. 105, 2002, 
pp. 539-542. 

[5] V. Chalana, D.T. Linker, D.R. Haynor, and Y. Kim, A multiple active 
contour model for cardiac boundary detection on echocardiographic se- 
quences, IEEE Trans. Med. Imaging (TMI), Vol. 15(3), 1996, pp. 290- 
298. 

[6] D. Comaniciu, X-S. Zhou, and S. Krishnan, Robust real-time myocardial 
border tracking for echocardiography: an information fusion approach. 
IEEE Trans. Med. Imaging (TMI), Vol. 23(7), 2004, pp. 849-860. 

[7] S. Goo, P. Joshi, G. Sand, D. Gerneke, A. Taberner, Q. DoUie, I. 

LeGrice, and D. Loiselle, Trabcculae carneae as models of the ventricu- 
lar walls: implications for the delivery of oxygen. Jour. Gen. Physiology, 
Vol. 134(4), Oct. 2009, pp. 339-350. 

[8] J. Kocnderink, Solid Shape, The MIT Press, Cambridge, Massachusetts, 
1990. 

[9] C. Li, C. Xu, C. Gui, and M.D. Fox, Level set evolution without re- 
initialization: a new variational formulation, Proc. IEEE Conf. CVPR, 
Vol. 1, 2005, pp. 430-436. 

[10] D. Lowe, Distinctive image features from scale-invariant keypoints, Intl. 
Jour. Computer Vision, Vol. 60(2), 2004, pp. 91110. 

[11] W. E. Lorensen and H. E. Cline, Marching cubes: A high resolution 3D 
surface construction algorithm. Computer Graphics, Vol. 21(4), 1987, 
pp. 163169. 

35 



[12] G.E. Mailloux, F. Langlois, P.Y. Simard and M. Bertrand, Restoration 
of the velocity field of the heart from two-dimensional echocardiograms, 
IEEE Trans. Med. Imaging, Vol. 8(2), 1989, pp. 143-153. 

[13] J. Matas, O. Chum, M. Urban, and T. Pajdla, Robust wide-baseline 
stereo from maximally stable extremal regions. Image and Vision Com- 
puUng Vol 22(10), 2004, pp. 761-767. 

[14] Medtronic Inc, The Visible Heart webpage: 

http: //www. visibleheart.com/index.shtml 

[15] F. Memoli and G. Sapiro, A theoretical and computational framework 
for isometry invariant recognition of point cloud data. Found. Comput. 
Math. Vol. 5(3), 2005, pp. 313-347. 

[16] M. Mignotte, J. Meunier and J.C Tardif, Endocardial boundary esti- 
mation and tracking in echocardiographic images using deformable tem- 
plates and Markov random fields. Pattern Analysis and Applications, 
Vol. 4(4), Nov. 2001, pp. 256-271. 

[17] N.J. Mitra, L.J. Guibas, J. Giesen, and M. Pauly. Probabilistic fin- 
gerprints for shapes, Proc. Symp. Geometry Processing, 2005, Vienna, 
Austria. 

[18] N.J. Mitra, L.J. Guibas, M. Pauly. Partial and approximate symmetry 
detection for 3D geometry, Proc. ACM SIGGRAPH, 2006, pp. 560-568. 

[19] MPEG-7 overview (version 10), ISO/IEC JTC1/SC29/WG11, Tech. 
Rep., 2004. 



36 



[20] A. Mukhopadhyay, Z. Qian, S.M. Bhandarkar, T. Liu and S. Voros, 
Shape analysis of the left ventricular endocardial surface and its ap- 
plication in detecting coronary artery disease, Proc. Intl. Conf. Fund. 
Imaging and Modeling of Heart (FIMH), 2011, New York, NY. 

[21] A. Mukhopadhyay, Z. Qian, S.M. Bhandarkar, T. Liu, S. Rienhart, and 
S. Voros, Morphological analysis of the left ventricular endocardial sur- 
face and its clinical implications, Proc. MICCAI , Nice France, Oct. 
2012. 

[22] A R. Osada, T. Funkhouser, B. Chazelle, and D. Dobkin, Shape distri- 
butions, ACM Trans, on Graphics, Vol. 21, No. 4, pp. 807-832, October 
2002. 

[23] M. Ovsjanikov, A.M. Bronstein, M.M. Bronstein and L.J. Guibas, Shape 
Google: A computer vision approach to invariant shape retrieval, Proc. 
NORDIA Workshop, 2009. 

[24] M. Ovsjanikov, Q. Mrigot, F. Mmoh, L.J. Guibas: One Point Isometric 
Matching with the Heat Kernel, Eurographics Symposium on Geometry 

Processing, 2010 

[25] Z. Qian, D. Metaxas, and L. Axel, Boosting and nonparametric based 
tracking of tagged MRI cardiac boundaries, Proc. MICCAI, LNCS 4190, 
2006, pp. 636-644. 

[26] J. Sivic, and A. Zisserman, Video Google: A text retrieval approach 
to object matching, Proc. IEEE Intl. Conf. Computer Vision (ICCV), 
2003, pp. 1-8. 

37 



[27] J. Sivic, and A. Zisserman Video Google: efficient visual search of videos, 
in Toward Category- Lev el Object Recognition, Springer LNCS, J. Ponce 
et al. (Eds.), Vol. 4170, pp. 127-144, 2006. 

[28] J. Sun, M. Ovsjanikov, and L.J. Guibas, A concise and provably infor- 
mative multi-scale signature based on heat diffusion, Proc. Symposium 
on Geometry Processing, 2009. 

[29] R. Toldo, U. CastelUani, and A. Pusiello. The bag of words approach for 
retrieval and categorization of 3D objects. The Visual Computer, Vol. 
26(10), Oct. 2010, pp. 1257-1268. 

[30] T. Zaharia, and F. Preteux, 3D Shape-based retrieval within the MPEG- 
7 framework, Proc. SPIE Conf. Nonlinear Image Proc. Pattern Anal. 
XII, Vol. 4304, 2001, pp. 133-145. 

[31] Y. Zhang and A.B. Hamza, Vertex-based diffusion for 3-D mesh denois- 
ing, IEEE Trans. Image Processing, Vol. 16(4), April 2007, pp. 1036- 
1045. 



38 



